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Abstract 

Phase ordering dynamics of the (2 + 1)- and (3 + l)-dimensional 
(j)^ theory with Hamihonian equations of motion is investigated nu- 
mericahy. Dynamic scahng is confirmed. The dynamic exponent z 
is different from that of the Ising model with dynamics of model A, 
while the exponent A is the same. 
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1 Introduction 



It is believed that macroscopic properties of many particle systems could 
be, in principle, described by microscopic deterministic equations of motion 
(e.g., Newton, Hamiltonian and Heisenberg equations), if all interactions, 
boundary conditions and initial states could be taken into account. However, 
practically it is very difficult to solve these equations, except for some simple 
cases. Therefore, statistical mechanics is developed to describe effectively 
statistical systems. In usual circumstances, ensemble theories give good de- 
scription of equilibrium states. But it is complicated for non-equilibrium 
states. A general theory does not exist. In many cases, some stochastic dy- 
namics, e.g., that described by Langevin-type equations of motion or Monte 
Carlo dynamics, is considered to be approximate theories. Anyway, whether 
microscopic equations of motion could really produce the results of statistical 
mechanics, or vice verse, remains open, e.g., see Refs H, ||, H, H. 

With recent development of computers, gradually it becomes possible to 
solve microscopic deterministic equations numerically. This attracts much 
attention of scientists in different fields. Study of microscopic fundamental 
dynamics is on the one hand to test statistical mechanics, and on the other 
hand to explore new physics. For example, assuming a system is isolated, 
there is only internal interaction, and periodic boundary conditions can be 
adopted. Computations are greatly simplified. To achieve ergodicity, the 
system should start from random initial states. Recently, such effort has 
been made for the 0{N) vector model and XY model 0]. The results 
support that deterministic Hamiltonian equations correctly describe second 
order phase transitions. The estimated static critical exponents are consis- 
tent with those calculated from canonical ensembles. More interestingly, the 
macroscopic short-time (non-equilibrium) dynamic behavior of the (2 + 1)- 
dimensional 0^ theory at criticality has also been investigated and dynamic 
scaling is found |^. The results indicate that Hamiltonian dynamics in 
two dimensions with random initial states is in a same universality class of 
Monte Carlo dynamics of model A. 

In a similar spirit, in Ref. [[l^] the author has investigated phase ordering 
dynamics of the (2 + l)-dimensional 0^ theory with Hamiltonian equations of 
motion. Assuming random initial states, there is a minimum energy density 
which is above the real minimum energy density of the system. Starting 
from this minimum energy density (noting that energy is conserved), which 
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is well below the critical energy density, phase ordering occurs. Dynamic 
scaling behavior is found. The dynamic exponent z is different from that of 
model A dynamics, but the exponent A governing the power law decay of 
the autocorrelation looks the same. Somewhat interesting is that the scaling 
function of the equal-time spatial correlation function is the same as that of 
the Ising model with model A dynamics. All above results are independent 
of the parameters in the system. 

The purpose of this paper is twofold: Firstly, we generalize the compu- 
tations to (3 + 1) dimensions. This is important since our realistic world is 
in (3 + 1) dimensions. Furthermore, in phase ordering of model A dynam- 
ics, the dynamic exponent z is dimension independent but the exponent A 
is dimension dependent. It is interesting whether this property is kept in 
Hamiltonian dynamics. Attention will be also put on whether the scaling 
function of the equal-time spatial correlation function in three dimensions 
is the same as that of the Ising model with model A dynamics. Secondly, 
to achieve more confidence on our conclusions we will reexamine the results 
for (2 + 1) dimensions in Ref. |Ty] using somewhat different, more careful 



approaches. Since the computations in (3 + 1) dimensions are very much 
time consuming, more accurate data are obtained in (2 + 1) dimensions. 

In the next section, we introduce the model and analyze the dynamic 
scaling behavior. In Sec. 3, numerical results are presented. Finally come 
the concluding remarks. 

2 Phase ordering dynamics 

In the following, we outline phase ordering dynamics with Hamiltonian equa- 
tions of motion. For a recent review of general ordering dynamics, readers 



are referred to Ref. [11 . 



2.1 The model 

For an isolated system, the Hamiltonian of the ((i + l)-dimensional (^^ theory 
on a square or cubic lattice is 



(1) 
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with TTi = (pi. It leads to the equations of motion 

4>i = + (Pi-f^ - 20*) + m^^t - -Lg^^i ■ (2) 

Here fi represents spatial directions. Energy is conserved in these equations. 
Solutions are assumed to generate a microc-anonical ensemble. The tempera- 
ture could be defined as the averaged kinetic energy. For the non- equilibrium 
dynamic system, however, total energy is a more convenient controlling pa- 
rameter, since it is conserved and can be taken as input from initial states. 
For given parameters and g, there exists a critical energy density ec, sep- 
arating the ordered phase (below ec) and disordered phase (above ec). The 
phase transition is of second order. 

We should emphasize that a Langevin equation at zero temperature is 
also 'deterministic' in the sense that there are no noises, but it is essentially 
different from the Hamiltonian equation (0). The former describes relaxation 
towards equilibrium at zero temperature for a non-isolated system, but the 
latter contains full physics at all temperatures for an isolated system. 

The order parameter of the 0^ theory is the magnetization. The time- 
dependent magnetization M = M^^\t) and its second moment M^"^^ are 
defined as 

r 1 (^) 

L is the lattice size and d is the spatial dimension. Here it is important that 
the average is over initial configurations. This is a real sample average and 
different from the time average in equilibrium. 



), A; = 1,2. (3) 



Following ordering dynamics with stochastic equations |11|, we consider 



a dynamic process that the system, initially in a disordered state but with 
an energy density well below ec, is suddenly released to evolve according to 
Eq. (^). For simplicity, we set initial kinetic energy to zero, i.e., 0j(O) = 0. To 
generate a random initial configuration {0j(O)}, we first fix the magnitude 
|0j(O)| = c, then randomly give the sign to 0j(O) with the restriction of a 
fixed magnetization in units of c, and finally the constant c is determined by 
the given energy. 

In case of stochastic dynamics, scaling behavior of phase ordering is dom- 
inated by the fixed point (T/,Tp) = (oo,0) with Tj being the initial tem- 



perature and Tp being the temperature after quenching |Tll|. In Hamilto- 



nian dynamics, the energy density can not be taken to the real minimum 
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Cmm = —3m^/2g since the system does not move. Actually, for the initial 
states described above, the energy is given by 



(4) 



For the case of d < it is demonstrated in Ref. |]T0[ that for a en- 

ergy density well below the critical point due to the competition of two 
ordered states, phase ordering occurs when the initial magnetization is set 
to zero. The scaling behavior is dominated by the minimum energy density 
Vmin = Vmin/L'^, which is a kind of fixed points. Above Vmin, there are extra 
corrections to scaling. From now, we redefine the energy density Cmm as zero. 
Then the fixed point is eo = fmm — Cmin- In this paper, we consider only the 
energy density at exactly the fixed point eo. 



2.2 Dynamic scaling behavior 

Let us first consider the case of the initial magnetization mo = 0. An impor- 
tant observable is the equal-time correlation function 

c{r,t) = ^^{j:<p.m+r{t)) . (5) 

The scaling hypothesis is that at the late stage of the time evolution, C (r, t) 
obeys a scaling form 

C{r,t) = f{r/t'/n , (6) 

where z is the so-called dynamic exponent. To the understanding of the au- 
thors, here 'late' is meant in microscopic sense. In other words, when the 
domain size (~ t^/') is big enough in units of the lattice spacing, scaling 
behavior emerges. At finite t, of course, there may be corrections to scaling. 
Corrections to scaling are generally not universal. They may induce difficul- 
ties for observing scaling behavior and uncertainties in the determination of 
the critical exponents. 

Simple understanding of the scaling behavior of C{r,t) can be achieved 
from the second moment of the magnetization. Integrating over r in Eq. (|^), 
we obtain a power law behavior 

M(2)(t) -t'^/^ . (7) 
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Another interesting observable is the auto-correlation function 

m = ^AT.<i^Mm)- (8) 

i 

The scaling hypothesis leads to a power law behavior 

A{t) ~ (9) 

which implies that ordering dynamics is in some sense 'critical'. Here A is 
another independent exponent. 

For the discussions above, the initial magnetization mo is set to zero. If 
mo is non-zero, the system reaches a unique ordered state within a finite time. 
If mo is infinitesimally small, however, the time for reaching the ordered state 
is also infinite and scaling behavior can still be expected, at least at relatively 
early times (in macroscopic sense). In this case, an interesting observable is 
the magnetization itself. It increases by a power law 

M{t)^t\ e = {d-x)/z. (10) 

The exponent 6 can be written as xo/z, with xo being the scaling dimension 
of mo. This power law behavior has deeply been investigated in critical 
dynamics |]T2|, 0. The interesting point here is that 6 is related to the 



exponent A which governs the power law decay of the auto-correlation. By 
combining measurements of 6 and A, one can also estimate the dynamic 
exponent z. 



3 Numerical results 

To solve the equations of motion (|]) numerically, we discretize (pi by {(j)i{t + 
At) + (piit — At) — 20j(t))/(At)^. Starting from an initial configuration, we 
update the equations of motion up to a certain maximum time t^ax- Then 
we repeat the procedure with other initial configurations. In Ref. |]TU[, rea- 
sonable results in two dimensions are obtained mainly with At = 0.05 up to 
tmax = 640 and a lattice size L = 521. 200 samples of initial configurations 
are used for averaging. For three dimensions, we also perform the computa- 
tions with At = 0.05 up to t^ax = 640 but with a lattice size L = 125. 50 
samples are taken for averaging. We have also carried out some computations 
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Figure 1: (a) C{r, t) in two dimensions obtained with L = 256 and At = 0.01 
is plotted with sohd hnes for t = 20, 40, 80, 160, 320, 640 and 1280 (from 
left). Circles fitted to the curve at the time t are the data at the time 2t but 
r being rescaled by a factor 2~^/^. C{r,t) obtained with At = 0.02 is also 
plotted with dashed lines in the figure but they overlap almost completely 
with the solid lines, (b) A{t) obtained with At = 0.01 (the solid line) in 
log-log scale. The curve for At = 0.02 overlaps completely with that for 
At = 0.01. 

with other At's and lattice sizes to confirm the results. At the time tmax, 
the equal-time correlation function C{r,t) decays to nearly zero at r ~ 45 
and this indicates also that the finite size effect with L = 128 is already 
small. Furthermore, in order to gain more confidence in our conclusions, es- 
pecially whether our t^ax has really reached the scaling regime, we perform 
more accurate computations in two dimensions (compared with those in Ref. 
0) with a lattice size L = 256 and At = 0.02, 0.01 up to tmax = 1280. 
The number of samples for averaging is 600. Somewhat different and more 
careful approaches will be adopted in this paper. 

In Fig. |l| (a), the equal-time correlation function C(r, t) in two dimensions 
is displayed. Solid lines are obtained with At = 0.01 and from left the time 
t is 20, 40, 80, 160, 320, 640 and 1280. Data for At = 0.02 are also plotted 
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Figure 2: a) C{r, t) in three dimensions obtained with L = 128 and At = 0.05 
is plotted with sohd hnes for t = 20, 40, 80, 160, 320 and 640 (from left). 
Circles fitted to the curve at the time t are the data at the time 2t but r being 
rescaled by a factor 2~^/^. The dashed line represents the scaling function in 
two dimensions, (b) A(t) and M^'^\t) in log-log scale. 



with dashed lines in the figures, but they almost completely overlap with the 
solid lines. For the curve oft = 1280, C{r, t) decays to nearly zero at r ~ 50. 
Therefore, we conclude that the finite size effect with the lattice size L = 256 
should be already negligible small. To confirm this, we have also compared 
the data with those in Ref . |T0[ . On the other hand, our data show that the 
finite At effect for At = 0.05 is also negligible. According to the scaling form 
(^, from data collapse of C{r, t) at different t's one can estimate the dynamic 



exponent z. As is observed in Ref. [10|, the effective dynamic exponent z{t) 



shows a small dependence on the time t. To explore this behavior and extract 
confidently the value of z, we perform scaling collapse of C{r, t) with the time 
t and 2t. In Fig. |l| (a), circles fitted to a solid line of the time t are the data of 
the time 2t with r being rescaled by a factor 2~-^/^, i.e., C(r, t) = C(r2^/^, 2t). 
The dynamic exponent z{t) is determined by the best fitting of the circles 
to the corresponding solid line. From the figure, we see clearly that the data 
collapse very nicely. 
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Figure 3: (a) The effective dynamic exponent z{t) measured from scaling 
collapse of C{r,t) with the times t and 2t. (b) Taking a = X/z, A{t)t'^ tends 
to a constant. 

In Fig. (a), a similar figure as Fig. ^ (a) is shown for C{r,t) in three 
dimensions. Scaling collapse is also observed, even though for larger r it is 
not as good as in two dimensions. This can be neither a finite size effect nor a 
finite At effect, since it exists also for small t's. To see the trend of z{t) as the 
time t evolves, in Fig. | (a) we plot the effective exponent z(t) against 1/t. 
For two dimensions, z{t) starting from a value around 3 gradually decreases 
and reaches 2.63(2) at t = 640 (i.e., obtained with data of C(r, t) at the time 
t = 640 and 2t = 1280). Assuming the behavior of z{t) will not changed 
essentially after t = 1280, the extrapolated value of z to the infinite time t 
is estimated to be 2.6(1). Interestingly, for three dimensions the exponent 
z{t) starting from a value around 2.5 increases slowly, but stabilizes at 2.7 
after t = 80. A good estimate of z is z = 2.7(1). Within statistical errors, 
the values of the dynamic exponent z in two and three dimensions coincide 
with each other, thus indicating that the dynamic exponent z is dimension 
independent. This can also be seen from the joining of two different curves 
at relatively larger times in Fig. ^ (a). 

In the case of the Ising model with Monte Carlo dynamics, the effective 
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exponent z{t) in two dimensions converges to z = 2 rather fast, e.g., see Ref. 



IJ] , but relatively slowly in three dimensions due to corrections to scaling. It 
might be somewhat general that phase ordering dynamics in three dimensions 
is somewhat more complicated than in two dimensions. 



An interesting fact observed in Ref. ||T0[ is that even though the dynamic 
exponent z of the 0^ theory in two dimensions with Hamiltonian dynamics is 
different from that of the Ising model with Monte Carlo dynamics, the scaling 
function f{x) in Eq. is the same. However, this is probably only by chance 
since it is not the case in three dimensions. The scaling function f{x) of the 
three dimensional 0"^ theory with Hamiltonian dynamics is different not only 
from that of the two dimensional but also from that of the three dimensional 
Ising model with Monte Carlo dynamics. The dashed line in Fig. ^ (a) shows 
the f{x) of the two-dimensional (p'^ theory. In general, Hamiltonian dynamics 
for isolated systems differs indeed from stochastic dynamics for non-isolated 
systems. 

For simple understanding of the correlation function C(r, t), one can mea- 
sure the time-dependent second moment M^^^ (t) . The scaling form results in 
a power law behavior for M^"^^ (t) and from the slope in log-log scale one can 
estimate the corresponding exponent. Such an approach is rather typical and 
useful in critical dynamics [|l^. It can also be applied in ordering dynam- 



ics, but less efficient. In critical dynamics, in the scaling collapse of C{r,t) 
one has to determine two exponents, the dynamic exponent z and the static 
exponent 2(3/ v. Therefore, it is efficient to read out directly the exponent 
{d — 2(3 / u) / z from the slope of of M'^'^\t) in log- log scale [1^. However, in 
ordering dynamics the 'static' exponent 2/?/z/ = and the scaling collapse of 
C(r, t) is only a one parameter fit. Measurements of M*^^^(t) do not show ad- 
vantage since it is not self-averaged and there is larger fluctuation for bigger 
lattices. This is seen from the data in Ref. |T^. Anyway, in Fig. |^ (b) we 
have plotted the second moment in log-log scale for the three dimensional (f)^ 
theory. The power law behavior is seen after t ~ 80 and this is consistent 
with Fig. (a). According to Eq. (|^), the resulting dynamic exponent is 
z = 2.5(2), consistent within errors with z = 2.7(1) measured from C{r,t). 

Another interesting exponent in ordering dynamics is A governing the 
power law decay of the auto-correlation A{t) in Eq. (H). The measurements 
of the auto-correlation in ordering dynamics is easier than in critical dynamics 
since the fluctuation is much smaller. The results for the 0^ theory in two 
and three dimensions are shown in Fig. |1] (b) and (b). In order to see 
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how the effective exponent Xj z depends on the time t, we have measured the 
slope of the curves in a time interval [t,2t]. The results are given in Table 
|I|. For both two and three dimensions, the exponent Xj z becomes stable 
after t = 160. The final values are X/z = 0.466(3) and 0.618(4) for two and 
three dimensions respectively. To show clearly that our estimates of X/z are 
indeed reasonable, in Fig. ^ (b) we plot A{t)t"' as a function of the time t. A 
correct value a = X/z should result in a constant for y4(t)t", at least for larger 
times. Such a behavior is nicely seen from the lower solid line and the dashed 
line for two and three dimensions in the figure. To confirm that the value 
X/z = 0.466(3) for two dimensions is really different from X/z = 0.625 for 
stochastic dynamics, the corresponding curve with a = 0.625 is also displayed 
there (the upper solid line). Obviously, it does not tend to a constant. 

t 40 80 160 320 640 

~Yd 0.508(1) 0.492(1) 0.469(7) 0.461(6) 0.463(6)~ 
3d 0.633(4) 0.609(1) 0.617(3) 0.619(7) 

Table 1: The exponent X/z measured in a time interval [t, 2t] from the auto- 
correlation in two and three dimensions. 

From measurements of z (from C(r, t)) and X/z, we estimate the exponent 
A = 1.21(5) and 1.67(6) for two and three dimensions respectively. For 
stochastic dynamics, theoretical prediction for two dimensions is A = 1.25 
[p!5| , PH , but in Monte Carlo simulations it is usually slightly smaller . 
Extrapolation is needed to obtain a value very close to 1.25. There is always 
some uncertainty in extrapolation. Therefore, we tend to claim that A of the 
(f)'^ theory in two dimensions with Hamiltonian dynamics is the same as that 
of stochastic dynamics. In three dimensions, our A = 1.67(6) agrees very well 



with the 'best' theoretical prediction 1.67 for stochastic dynamics [p!6| , |ll 
Numerical measurements of A for stochastic dynamics in three dimensions 
look somewhat problematic and the results fluctuate around the theoretical 
values. 

To complete our investigation, we have also simulated the initial increase 
of the magnetization in Eq. (plQl). Since the exponent 6 is relatively big, com- 
pared with that in critical dynamics |10, 13], we need to prepare a very small 



initial magnetization mo. In Fig. |^, the magnetization in three dimensions is 
plotted in log-log scale for mo = 0.00123, 0.00245 and 0.00491 (from below) 
respectively. The power law behavior is observed after t ~ 50. From the 
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slope, we measure the exponent 6. Within statistical errors, we can not find 
any mo dependence of 9. The value of 9 is estimated to be 0.55(2). With 9 
in hand, combining X/z = 0.618(4), we obtain another value for the dynamic 
exponent, z = 2.6(1). 

In Table U, all the exponents measured for the 0^ theory with Hamilto- 
nian dynamics are summarized. Results for two dimensions are taken from 
Ref. [0, but X/z, X and z from C{r,t) are slightly modified according to 
new data in this paper. Different measurements in two and three dimen- 
sions suggest that z = 2.6(1) is a good estimate for the dynamic exponent. 
Different from the case of critical dynamics |]13[, the critical exponent 9 in 



phase ordering dynamics has not yet drawn enough attention, even though 
it has been addressed [^. One reason might be that in ordering dynamics, 
increasing of the magnetization is expected if a non-zero initial value mo is 
set, but in critical dynamics, this is anomalous. Anyway, we think 9 is inter- 
esting since it gives another independent estimate for the dynamic exponent 
z or A. 







z 






9 X/z 


d/{X/z + 9) C{r,t) M^''^ 


A 


2d 


0.31(1) 0.466(3) 


2.6(1) 2.6(1) 2.6(1) 


1.21(5) 


3d 


0.55(2) 0.618(4) 


2.6(1) 2.7(1) 2.5(2) 


1.67(6) 



Table 2: Exponents of the 0^ theory with Hamiltonian dynamics. To calcu- 
late A, z measured from C{r,t) is taken as input. 



4 Conclusions 

In conclusion, we have numerically solved the Hamiltonian equations of mo- 
tion for the two- and three-dimensional 0^ theory with random initial states. 
Phase ordering dynamics is carefully investigated. Scaling behavior is con- 
firmed. The dynamic exponent z is dimension independent. Different mea- 
surements yield a value z = 2.6(1) and it is different from z = 2 for stochastic 
dynamics of model A. The scaling function for the equal-time spatial correla- 
tion function is dimension dependent, and in general also different from that 
of stochastic dynamics of model A (it is the same probably only by chance in 
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two dimensions). However, the exponent A of Hamiltonian dynamics is the 
same as that of stochastic dynamics of model A. 
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Figure 4: The magnetization in three dimensions in log-log scale. The lattice 
size is L = 128. Prom below, mo = 0.00123, 0.00245 and 0.00491. 
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